Magnetic properties of lightly doped antiferromagnetic YBa2Cu30j/. 
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The present work addresses YBa2Cu3 0j; at doping below x — 6% where the compound is a 
collinear antiferromagnet. In this region YBa2Cu30y is a normal conductor with a finite resistivity 
at zero temperature. The value of the staggered magnetization at zero temperature is « O.GfiB, 
the maximum value allowed by spin quantum fiuctuations. The staggered magnetization is almost 
independent of doping. On the other hand, the Neel temperature decays very quickly from Tjv = 
42074' at X = to practically zero at a; ~ 6%. The present paper explains these remarkable properties 
and demonstrates that the properties result from the physics of a lightly doped Mott insulator with 
small hole pockets. Nuclear quadrupole resonance data are also discussed. The data shed light on 
mechanisms of stability of the antiferromagnetic order at a; < 6%. 

PACS numbers: 74.72.-h 75.25.+Z 76.75. +i 78.70.Nx 
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I. INTRODUCTION 

It is well known that cupratcs arc layered compounds 
consisting of Cu02 planes and there are no doubts that 
the generic physics of cuprates are related to the Cu02 
plane. In spite of the same generic physics specific prop- 
erties of cuprates can be very different depending on crys- 
tal structure, ways of doping etc. The goal of the present 
work it to shed light on the generic physics via under- 
standing of specific properties of lightly doped antiferro- 
magnetic YBasCugOj; (YBCO). 

Cuprates are essentially doped Mott insulators. It is 
well established that a Mott insulator possesses a long 
range antiferromagnetic (AF) order, therefore, one of the 
generic problems is how the AF order evolves with dop- 
ing. Another generic problem is the shape of the Fermi 
surface. Are there small hole pockets as one expects 
for a very lightly doped Mott insulator, how the surface 
evolves with doping? 

Cuprates are intrinsically disordered materials because 
of mechanisms of doping. Disorder complicates a theo- 
retical analysis of experimental data usually masking the 
generic physics. YBCO is probably the least disordered 
cuprate in the low doping regime. In this paper I denote 
the hole concentration per unit cell of the Cu02 layer by 
X, this is the "doping". YBCO is not superconducting 
below X ~ 0.06 where it behaves as a normal conductor 
with delocalized holes. The zero temperature resistivity 
remains finite^ apart of a very weak logarithmic temper- 
ature dependence^i^ expected for a weak disorder. The 
heat conductivity also indicates delocalization of holesi^ 
This is very much different from La2-2;Srj;Cu04 where 
holes are localized and hence the compound is an Ander- 
son insulatori^ at x < O.lSfi Ultimately, at the very low 
doping, X < 0.01, the disorder wins even in YBCO and 
it also becomes the Anderson insulatorii It is helpful to 
have in mind an approximate empiric formulaii^ 
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The static "staggered" magnetization in YBCO has 
been recently measured in the ^SR experiment!^ The ex- 
perimental plot of the zero temperature magnetization 
versus doping is shown in Fig. [TJtop). Remarkably the 
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to relate the doping level x and the oxygen content y in 
underdoped YBCO at x < 0.12. 



FIG. 1: (Color online). The YBCO ^iSR data from 7?e/. Q. 
Top; Zero temperature "staggered" magnetization versus dop- 
ing. The solid line shows results of the present calculations. 
Bottom; The Neel temperature and the staggered magneti- 
zation versus doping. The inset shows the superconducting 
critical temperature versus doping. Large red dots show the 
Neel temperature calculated in the present work. 

zero temperature magnetization is almost doping inde- 
pendent up to a: sa 0.06 and then it quickly decays. It 
is known from neutron scattering experiment»i2r— that 
the static magnetization fully disappears at the Quan- 
tum Critical Point (QCP) x « 0.09 indicating transition 
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to a state without static magnetism. Importantly, the 
magnetism at x > 0.06 is incommensurate, this is why 
in the first sentence of this paragraph I put "staggered" 
in inverted commas. Value of the incommensurate wave 
vector Q divided by 27r versus doping is plotted in Fig. [2] 

0.15| 
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FIG. 2: (Color online). Incommensurate wave vector versus 
doping. The blue square^, the red circle^, and the red dia- 
monS^ show neutron scattering data. The solid line shows 
the theoretical valuei^ 

While in the collinear antifcrromagnctic phase the zero 
temperature staggered magnetization is almost indepen- 
dent of doping, the Neel temperature decays very quickly 
from Tn = 420-ftr at x = to practically zero at a; « 6%. 
This is shown in Fig. [TJbottom) copied from Ref. H. The 
present paper explains these puzzling magnetic proper- 
ties and shows that they are related to small hole pockets 
of lightly doped Mott insulator. 

YBa2Cu30y is doped via filling oxygen chains located 
above the Cu02 planes. It has been argued that at 
y = 6.5, where the every second chain is full, the chain 
modulation causes the corresponding charge density wave 
(CDW) of in-plane holesJ^ii^ Nuclear quadruple reso- 
nance (NQR) for in-plane Cu is an excellent local probe 
of the hole dcnsityJi Fig 13] shows ^'^Cu NQR frequency 
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FIG. 3: ( Color online). In-plane Cu NQR frequency sweeps 
from Ref. [73. 

sweeps from Ref. [l^ for several values of oxygen con- 
tent. There is a single narrow line at about optimal 
doping y ~ 1 indicating a very homogeneous hole den- 
sity corresponding to completely filled chains. On the 
other hand, at y « 6.5 there are two distinct lines indi- 
cating a bimodal hole density distribution in agreement 
with Refs. [T5IIT6I . Importantly, the bimodal distribution 
is evident even at lower doping, y = 6, 4, 6.45, indicating 



the CDW induced by the oxygen chain superstructure. 
Below y = 6.5 the NQR lines are broader compared to 
y = 6.5. This is because the oxygen superstructure with 
every second chain filled cannot be perfect away from 
y = 6.5. It is worth noting that holes go to the Cu02 
plane only at y > 6.2, sec Eq.([T]). In the undoped case, 
X = Q, y < 6.2, there is only one NQR line with frequency 
vq w 23.3MHz independent of y, sec Ref. [I^ Compar- 
ing vq with frequencies of lines in FiglS] we see that the 
hole doping shifts the NQR frequency very stronglyJi 
The present paper explains significance of the chain in- 
duced CDW for stability of the collinear antiferromag- 
netic phase at a; < 0.06. 

Structure of the paper is the following. The effective 
theory describing YBCO at low doping was formulated 
before in Ref. Il4l . Section II summarizes ideas of the 
effective theory. In Section III the theory is applied to 
calculate reduction of the staggered magnetization in the 
antifcrromagnctic phase at zero temperature. Temper- 
ature reduction of the staggered magnetization at zero 
and nonzero doping, < x < 0.06, is calculated in Sec- 
tion IV. Interplay between the chain induced CDW, small 
hole pockets, and stability of the collinear antiferromag- 
netic phase is discussed in Section V. Section VI presents 
conclusions of the paper. 



II. EFFECTIVE LOW ENERGY THEORY 
DESCRIBING LIGHTLY DOPED YBCO 



This section summarizes the most important points of 
the effective low energy suggested in Refs il4li20l to de- 
scribe YBCO at low doping. The analysis is based on 
the two-dimensional t — t' — t" — J model at small dop- 
ing. The generic case of the single layer has been con- 
sidered in Ref. [2^ After integrating out the high energy 
fluctuations one comes to the effective low energy ac- 
tion of the model. The effective low-energy Lagrangian 
is written in terms of the bosonic n-field (n^ = 1) that 
describes the staggered component of the copper spins, 
and in terms of fermionic holons tp. The term "holon" 
is used instead of "hole" to stress that spin and charge 
are to large extent separated, see Ref. [20. The holon 
has a pseudospin that originates from two sublattices, 
so the fermionic field tJj is the spinor in the pseudospin 
space. Minimums of the holon dispersion are at the nodal 
points qo = (±7r/2, ±7r/2). So, there are holons of two 
types corresponding to two pockets. The dispersion in 
a pocket is somewhat anisotropic, but for simplicity let 
us use here the isotropic approximation, e (p) sa ^/3p^ , 
where p = q — qo- The lattice spacing is set to be equal 
to unity, 3.81 A— > 1. All in all, the effective Lagrangian 
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for the single layer reads2£ 



y(Vn) 



(2) 



The first two terms in the Lagrangian represent the usual 
nonlinear cr model. The magnetic susceptibility and the 
spin stiffness are x± ~ 0.53/8 « 0.066 and ps ~ 0.175^. 
Hereafter the antiferromagnetic exchange of the initial 
t-J model is set to be equal to unity, 



J w 130 meV 



1. 



Note that ps is the bare spin stiffness, therefore by def- 
inition it is independent of doping. The rest of the La- 
grangian in Eq. ([2]) represents the fermionic holon field 
and its interaction with the n-field. The index a — a^b 
indicates the pocket in which the holon resides. The pseu- 
dospin operator is ^a, and Gq — {±1/^/2, l/\/2) is a unit 
vector orthogonal to the face of the magnetic Brillouine 
zone (MBZ), where the holon is located. The argument 
of in Eq. ([2]) and the time derivative of the fermionic 
field in the same equation are "long" (covariant) deriva- 
tives, 

V = -iV + i<T • [n X Vn] 

T^t = dt + -(7 ■ [n X ft] . 

The covariant derivatives refiect gauge invariance of the 
initial t — t' — t" — J model. 

Numerical calculations within the t — t' — t" — J model 
with physical values of hopping matrix elements give the 
following values of the coupling constant and the inverse 
mass, g ~ 1, /3 « 2.4. The value of the inverse mass 
/? = 2.4 corresponds to the effective mass m* = l.Sme- 
The dimensionless parameter 



A: 



(3) 



plays the defining role in the theory. If A < 1 , the ground 
state corresponding to the Lagrangian ^ is the usual 
Neel state, the state is coUinear at any small doping. If 
1 < A < 2, the Neel state is unstable at arbitrarily small 
doping and the ground state is a static or a dynamic spin 
spiral. The wave vector of the spiral is 



Q = —X 

Ps 



(4) 



If A > 2, the system is unstable with respect to phase 
separation and/or charge-density- wave formation and 
hence the effective long-wave-length Lagrangian ([2|) be- 
comes meaningless. The pure t — J model [t' = t" = 0) 
is unstable since it corresponds to A > 2. 



To find parameters of the effective action ([2]) one can 
rely on calculations within the t — t' — t" — J model or 
alternatively one can fit experimental data. Both ap- 
proaches produce very close values of the parameters. 
The fit of elastic and inelastic neutron scattering data for 
La2-2;Sr2;Cu04 performed in Ref. [lO gives the following 
values, 5 = 1, /3 sa 2.7 (m* = 1.5me), A « 1.30. The fit 
of data on magnetic quantum oscillations in YBa2Cu30y 
performed in Ref. [l^ gives two possible sets, 



g = l, 13 = 2.78 (m* 
c, = 1, 13 = 2.95 (to* 



1.45TOe), A =1.31, 
1.35TOe), A = 1.23 . (5) 



These values will be used in the present work. 

It is very easy to understand the reason for instability 
of the commensurate AF order under doping. Assum- 
ing such an order one can calculate the magnon Green's 
function 



G{uj,q) cx 
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uj^ — c^q^ — V{uj, q) + iO 



(6) 



where c = \/ Ps/X-L ~ 1.17-\/2J is the magnon speed 
in the parent Mott insulator and V{uj,q) is fermionic 
polarization operator. A well known peculiarity of the 
two-dimensional (2D) polarization operator is its inde- 
pendence of doping as soon as w = and q is sufficiently 
small. A straightforward calculation gives at g 0, 
V{0,q) = —Xc^q^. Hence, at A > 1 the Stoner criterion 
in ([6]) is violated, the Green's function possesses poles 
at imaginary frequency indicating instability of the AF 
ground state at an arbitrary small doping. 

In YBCO the AF order is commensurate at a; < 0.06, 
therefore the effective action ([2]) cannot be directly ap- 
plied to this compound. To understand YBCO one can 
certainly assume that A is doping dependent, A < 1 at 
X < 0.06 and A > 1 at a; > 0.06. Purely theoretically 
it is hardly possible to have a significant cc-dependence 
of A, but as a scenario one can consider this. However, 
in this scenario the incommensurate wave vector Q must 
jump from Q = at a; < 0.06 to Q given by Eq.Q at 
X > 0.06. This is not consistent with data, there is no a 
jump, the incommensurate wave vector evolves smoothly 
above x = 0.06, see Figl^l 

A model describing the smooth evolution of Q with 
doping was suggested in Ref. [13. In addition to ([2]) 
the model incorporates two points, (i) Due to the bi- 
layer structure the magnon spectrum in YBCO is split 
into acoustic and optic mode. The optic gap is about 
70meVr^ hence the optic mode does not influence the 
low energy dynamics, only acoustic magnons are impor- 
tant for these dynamics, (ii) The second point of the 
model is an assumption that the fermionic dispersion is 
split in two brunches as it is shown in Fig Hljleft). The 
splitting is Aq. The effective action that originates from 
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FIG. 4: (Color online). Filling of split holon hands m YBCO 
at X < xo (left) and x > xo (right). The solid and the dashed 
lines in the spin-spiral state, x > xo, correspond to different 
pseudospin projections, the splitting is igQ. The doping xi 
indicated on the top of the right figure is xi = 0.5xo /(A — 1), 
see Ref UI- 



(21) and incorporates these two points reads 



ct—a,b 7 — ±1 



•r a,' 



ea(^) -7— 



(7) 



Compared to ^ the first line is multiphcd by two since 
the bilayer has the twice larger spin stiffness and mag- 
netic susceptibility. In addition to the pocket index a 
the holon field 'il'a,-y gets an additional index 7 ~ ±1 
that indicates the branch of the split dispersion as it is 
shown in Fig U) Originally the paper— suggested that 
the hole band splitting 7 = ±1 was due to the hole hop- 
ping between layers inside the bilayer. So, Aq was the 
bonding-antibonding splitting. However, our recent anal- 
ysis^i indicates that antiferromagnetic correlations forbid 
the bonding-antibonding splitting. So, contrary to the 
assumption in Ref. [13 the interlayer hopping cannot con- 
tribute to Aq. In Section VI of the present paper I argue 
that the splitting Aq is due to oxygen chains. For now 
let us accept the action ([T]) and study consequences of 
this action. 

When doping is sufficiently small, 



X < Xq 



Ao 

7r/3 ' 



(8) 



only the 7 = 1 band is populated, see Fig. [Deleft). In 
this case the fermionic polarization operator is a half of 
that for the single layer case, 7^(0, q) = —^Xc^q'^, g — > 0. 
Hence the Stoner stability criterion in Eq. (jG)) is fulfilled 



and the Neel order is stable. According to both neutron 
scattering datapiSiiiiH shown in Fig. [5] and to /iSR datai^ 
shown in Fig. [T] the value of xq is xq sa 0.06. Hence, due 
to Eq.® the band splitting is 

Ao « 0.5J« 65meF . (9) 
At X > Xo fermions populate both 7=1 and 7 = — 1 
bands. Fig. fright), the polarization operator is doubled 
compared to the x < xq case, and the Stoner instability 
is there. As a result at x > xq the system develops the 
spiral with the wave vector— 



g X -Xq 

Ps3-2X 



(10) 



The plot of Q/27r versus doping is shown in Fig. [2]by the 
solid line. The development of the spin spiral is driven by 
the pseudospin splitting of the fermionic bands ±gQ as it 
is shown by solid and dashed lines in Fig. fright). Thus, 
Xq is a Lifshitz point, where the 7 = — 1 band starts to 
populate, and where simultaneously the spin spiral starts 
to develop. In the present paper I consider quantum and 
thermal fluctuations in the Neel state, x < xq. Quantum 
fluctuations in the spin spiral state sd, x > xq will be 
considered separately<^ 

To summarize this section. Small hole pockets and 
associated spin spiral state are generic properties of all 
cuprates at low doping. The key point of the YBCO phe- 
nomenology is splitting of the hole pockets. This splitting 
together with splitting of magnon to the acoustic and the 
optic mode provides stability of the AF order up to 6% 
doping. 



III. QUANTUM FLUCTUATIONS IN THE 
NEEL PHASE, REDUCTION OF THE ZERO 
TEMPERATURE STAGGERED 
MAGNETIZATION. 

There are two magnons (two polarizations) in the Neel 
phase a,t X < Xq. The Green's function of each magnon 
reads 



G(c.,q) = 



(2X±)-^ 



cV-^o(w,q) + *0 



(11) 



Only the 7 = 1 band, see Fig. IHleft), contributes to 
the polarization operator Voicdjq). Calculation of the 
polarization operator for the single layer was performed 
in Ref. [l^. Comparing the single layer action ([2]) and 
the double layer action ([7]) and having in mind that at 
X < Xq only the 7 = 1 band is occupied one immediately 
concludes that in the double layer case the polarization 
operator is a half of that calculated in Ref. [20. Hence 
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V 



(12) 



Here -pp is the Fermi momentum of the 7 = 1 band and 
6{x) is the usual step function. I've aheady pointed out 
above that at q < 2pp, Vo{0,q) = —^c^q"^, so the Neel 
state is stable if A < 2. 

It is instructive to look at the magnon spectral function 
—ImG{LO,q) that describes inelastic neutron scattering. 
Spectral functions for x = 0.05 and for three values of the 
momentum q arc plotted in Fig.[5]by solid lines. Spectral 
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x=0.05, pF=0.56 
q/2;t=0.1*pF/27l =0.009 



q/27i=0.2*pF/27i =0.018 

q/2ji=0.4*pF/27i=0.036 
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FIG. 5: (Color online). Solid lines show magnon spectral 
functions —ImG{uj,q) versus w for three values of momen- 
tum q and for doping x = 0.05. Dashed lines show spectral 
functions m the parent Mott insulator at the same values of 
momentum. 

functions for both sets of parameters from Eq. ([S]) are al- 
most identical. To be specific I present functions for the 
second set. In the same Fig.[5]thc dashed lines show spec- 
tral functions in the parent Mott insulator [Po(w, q) — 0]. 
The spectra demonstrate the low energy incoherent part 
that absorbs more that 50% of the spectral weight. The 
magnon quasiparticle peaks are still clearly pronounced. 
Their intensities are about half of that in the parent com- 
pound, and positions arc slightly shifted up compared 
to the parent, the shift is proportional to the doping, 
5(jjq oc X. It is worth noting that while the reduction of 
the spectral weight is a reliable result, the upward shift is 
probably a byproduct of the low energy effective theory. 
The effective theory accurately accounts for the magnon 
"repulsion" from the particle-hole continuum that is be- 
low the magnon. The "repulsion" results in the upward 
shift. However, there is also a "repulsion" from very high 
energy excitations {E ^ 2t ^ 6 J) that are related to the 
incoherent part of the hole Green's function. This re- 



pulsion, unaccounted in the effective theory, leads to the 
downward shift of the magnon frequency that is also pro- 
portional to doping)^ More generally one can say that 
the chiral effective field theory employed in the present 
work allows controllable calculations of effects that are 
X- independent or scale as yG; or a;ln(x). Quantities that 
scale as the first or as a higher than first power of x are 
generally beyond the scope of the theory. Therefore, at 
this stage one can say only that position of the magnon is 
approximately the same as that in the parent compound, 
but the magnon spectral weight is significantly reduced. 
Another point worth noting is the absence of the hour- 
glass dispersion. The low energy incoherent part of the 
Green's function clearly pronounced in Fig. [5] is trans- 
formed to the coherent hourglass only at x > xq, beyond 
the Lifshitz point 1^ 

Quantum fluctuation of the staggered magnetization is 
given by the standard formula 

K)--2E / ^G(.,q)=-2^ / ^/™G(.,q) 
q q 

(13) 

The factor 2 comes from two polarizations. This ex- 
pression must be renormalized by subtraction of the 
ultraviolet-divergent contribution that corresponds to 
the undoped c-modcl. The integral in (|T3|) can be calcu- 
lated analytically with logarithmic accuracy 
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(14) 



There are two points to note, (i) In spite of the ultravio- 
let renormalization (cr-model subtraction) the fluctuation 
depends on the ultraviolet momentum cutoff A ~ 1. (ii) 
The leading logarithmic term, a;ln(A^/a;), comes from 
momenta pi? <C g <C A. 

The logarithm ln(A^/a;) is not large, the logarithmic 
accuracy is not sufficient. Fortunately a numerical in- 
tegration of (|13p is straightforward. The result is pre- 
sented in Fig. |6l where (n^) is plotted versus doping. 
The second set of parameters from jl]) is used, results 
are presented for two values of the ultraviolet cutoff 
A. Reduc tion of the static component of the n-field is 

i(ri^). Hence the staggered 
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magnetization is 
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(15) 
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FIG. 6: (Color online). Quantum fluctuation {n\) versus 
doping for two values of the ultraviolet cutoff A. 



Here I take into account that the used regularization 
procedure corresponds to the normaUzation of the static 
component of the n-ficld to unity at zero doping when 
the staggered magnetization is 0.615/XB<2i The plot of 
the calculated staggered magnetization fi versus dop- 
ing together with experimental data^ is presented in 
Fig. [Ijtop). Dependence of the theoretical curve on A 
is pretty weak, to be specific the curve corresponding to 
A = 1 is presented. Agreement between the theory and 
the experiment in the Neel phase is excellent. Thus, it 
is understood why quantum fluctuations only slightly re- 
duce the staggered magnetization. 

Note that the presented calculation is valid only in the 
Neel phase, x < 0.06. Physics in the spin-spiral phase, 
X > 0.06, is very much different because of the appear- 
ance of the soft "hourglass" dispersion and consequently 
because of greatly enhanced quantum fluctuations. The 
corresponding results will be published scparatelyi^ 



IV. TEMPERATURE DEPENDENCE OF THE 
STAGGERED MAGNETIZATION IN THE NEEL 
PHASE 

A. Zero doping 

It is well known that due to the Mermin- Wagner theo- 
rem the Neel temperature in a spin-rotationally-invariant 
2D system is exactly zero, Tjv = 0. Cuprates are lay- 
ered systems with a very small Heisenberg interaction, 
J_L < 10"'' J, between layers or bilayers. In spite of 
its smallness the interaction makes the system three di- 
mensional and hence it makes the Neel temperature fi- 
nite, Tjy ^ J/ \n{J/ J±). Temperature dependence of the 
staggered magnetization in layered Heisenberg antiferro- 
magnets has been intensively studied theoretically, for a 
review see Ref. [l^- Unfortunately there is no a "small 
theoretical parameter" in the problem, therefore, while 
a qualitative behaviour is absolutely clear, there is no 
a universal quantitative description, different theoretical 
approaches give quite different results.— In the present 
section I develop an effective description of the tempera- 
ture dependence. This description certainly is not a rig- 
orous solution of the layered Heisenberg antiferromagnet 



for all temperatures. This is a sort of interpolation be- 
tween T <g; Tn regime and T « Tjv regime. Importantly, 
the "interpolation" allows to describe quantitatively an 
undopcd layered Mott insulator, and much more impor- 
tantly it allows to move to the finite doping in the next 
subsection. 

Let us start from the single layer case (La2Cu04) and 
rewrite Eq. p^ in the Matsubara technique at a finite 
temperature. 



(ni) 
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X- 
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(16) 



where 



cq and = 2sttT, s = 0,±1,±2,... is the 



Matsubara frequency. Hence equation for — 1 — ^ (n-'j) 
can be rewritten in the rcnormalization group (RG) form 



T 



nzdln{q) 27: p. 



(17) 



where psq — Psiq) is the g-dcpcndcnt spin stiffness. 
Eq. ()17|) assumes 2D geometry, so it is valid at q > qmin, 
where the infrared cutoff qmin VJT is due to the 
Heisenberg coupling along the third dimension. To solve 
the RG problem one needs to add information how the 
spin stiffness scales with the magnetization. Let us write 
the relation between the magnetization and the spin stiff- 
ness as 



dp, 



sq 



(18) 



It is known^ that one loop calculation valid at ti^ « 1 
results in r = 1 that implies p (x Uz- On the other 
hand, close to the Neel temperature when <C 1 one 
should expect scaling very close to quadratic, p oc (r « 
2). This is because the critical index 77 of the magnon 
quasiparticle residue is very small, see, e.g. Rcfs. |28ll29|. 
For now I keep the power r as a parameter. Eqs. (jl7|) 
and psp are combined to 



dp. 



sq 



dln{q) 
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(19) 



To perform the ultraviolet rcnormalization let us intro- 
duce PA, the spin stiffness at the ultraviolet cutoff. Then 
due to Eg. p^ the finite temperature spin stiffness at 
q = reads 



PsT = PA 



rT 
2^ 



E: 



■es 



(20) 



This expression can be renormalized by the condition 
that at zero temperature (more accurately at T <C cqmin) 
the spin stiffness is equal the standard value pso ~ 0.175 J 
corresponding the cr-modcl originated from the spin 1/2 
Heisenberg model. After the rcnormalization Eg. d^Ul) is 
transformed to 



PsT — PsO 



rT 
2^ 



^q 

2T 



dq 



(21) 
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FIG. 7: Staggered magnetization versus temperature in 
La2Cu04, left, and YBa2Cu306, right. In the left plot red 
circles show neutron scattering data.— In the right plot red 
circles show neutron scattering datai^ and green squares show 
/iSR data.— Theoretical curves with r = 2 are shown by solid 
lines and theoretical curves with r = 1 are shown by dashed 
lines. 



The 3D interaction J± fixes value of tlic infrared cutoff 
Qmin, however, one has to remember about scaling of the 
cutoff with the staggered magnetization n^, see Refill, 



QminO V n 



where, due to p^. 



PsT 
PsO 



l/r 



(22) 
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Eqs. (PT|) . (|22p . (l^51) can be easily integrated numeri- 
cally. The Neel temperature is determined by zero of 
the spin stiffness (HJ). The infrared cutoff q,ninQ is the 
only free parameter in the theory. Value of the pa- 
rameter has to be tuned up to reproduce the measured 
Neel temperature. It has to be clear that qmino orig- 
inates not only from J±, relativistic anisotropics like 
Dzyaloshinskii-Moria etc, also contribute to qminO- Let 
us recall that due to the used regularization procedure 
the staggered magnetization is /i = 0.615//b't-2, where 
0.615/xb is the the staggered magnetization in the par- 
ent Heisenberg modeli^ Staggered magnetization ver- 
sus temperature in La2Cu04 is presented in Fig. [Tljleft). 
Red circles show neutron scattering data<^ The theoret- 
ical curve with r = 2 and qminO = 0.024 is shown by 
the solid line and the theoretical curve with r = 1 and 
qminO = 0.004 is shown by the dashed line. The curve 
with r = 1 corresponding to the single loop RG describes 
the data very poorly. This illustrates the known problem 
of poor accuracy of the single loop RGi^ However, the 
curve with r = 2 corresponding to the critical scaling of 
the spin stiffness describes the data very well. 

In the double layer case (YBCO) the coefficient ^ 
before the integral in Eq. ipTj) has to be replaced by 
the twice smaller one, j^. The point is that the op- 
tic magnon in YBCO has a gap 70meV and therefore it 
does not contribute to the low energy dynamics. Only 
acoustic magnon is important, hence the effective num- 
ber of magnons is twice smaller compared to LCO. Neu- 
tron scattering datai^ for YBa2Cu306 are shown in 
FigEl^right) by red dots. Green squares show /iSR data^S 
The theoretical curve with r — 2 and qminO = 0.0085 is 



shown by the solid line and the theoretical curve with 
r = 1 and qminO = 0.0004 is shown by the dashed line. 
Again, the curve with r = 1 is not consistent with the 
data. The curve with r = 2 is quite good. 

It is worth noting that for both LCO and YBCO the 
values of the infrared cutoff qminO for r = 1 are unrealisti- 
cally small reflecting the same difficulty of single loop RG, 
see alsoj^ On the other hand, the cutoff values for r = 2 
are quite reasonable, they indicate that the Neel temper- 
ature is determined by spin-wave dynamics at distances 
up to 1/qminO ~ 100 lattice spacing along the plane. 

All in all, the conclusion is that the effective RG de- 
veloped in this subsection describes undoped compounds 
pretty well. To achieve this description one needs to set 
r = 2, this corresponds to the critical scaling of the spin 
stiffness expected in the vicinity of the Neel temperature, 
p (X n1. In the next subsection the developed description 
will be applied to the nonzero doping case. 



B. Nonzero doping 

To extend to the finite doping case one has to introduce 
in Eq. pip the fermionic polarization operator 




(24) 
dq 



Voiles, q) 2T q 



where the polarization operator Vq (i^^ , q) is calculated at 
Matsubara frequencies. Expression for the polarization 
operator follows from the Lagrangian ([T]). One can use 
vertexes derived in Ref. [l^ for the single layer case and 
rescale the vertexes by the factor 1 / ^/2 that follows from 
comparison of 1^ and (I7|). The polarization operator 
reads 



P - f 



p+q 



7=±1 p 



2T:\c^Pq^Re 



Jp 



(25) 



Here is the Fermi-Dirac distribution function 

1 



/■7 
Jp 



g(ep-7Ao/2-A')/T ^ ^ 



(26) 



with chemical potential /i (do not mix it up with mag- 
netic moment). Note that at T 7^ the 7 = — 1 band is 
also populated, see FigUJ^left). This is why the summa- 
tion in (|25p is performed over both bands, 7 = ±1. The 
chemical potential is determined by the condition 



2x=2x2y: e/; 

7=±1 p 



(27) 



that accounts for the double layer, the two pockets, and 
for the two psedospin projections. It is easy to check that 
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at zero temperature and at g < 2pp the zero frequency 
polarization operator is 7^0(01 9) = ~^^qj agreement 
with the real frequency analysis a.t x < xq in section III. 

Numerical evaluation of the polarization operator (|25p 
is straightforward. Substitution of the polarization oper- 
ator in the RG equation ([M)) and solution of this equa- 
tion together with ([^ and ([^ gives staggered mag- 
netization at a given doping and temperature. The RG 
equation is solved with r = 2 and qmino ~ 0.0085 as 
it has been discussed in the previous subsection. These 
parameters are relevant to the n-field and they are in- 
dependent of doping. Fermionic polarization operator 
([25|l is not very sensitive to the choice of parameters, 
to be specific I present results corresponding to the sec- 
ond set of parameters in Eq. ([S]). The band splitting 
Aq is determined by Eq.® that is responsible for the 
position of the Lifshitz point, xq — 0.06. Plots of the 
staggered magnetization versus temperature for several 
values of doping arc shown in Fig. [51 Theoretical curves 



0.6 















0.06 


iin'io 


14 0.03 \o.02 


0.01 








80 160 240 320 400 ^ 480 



FIG. 8: Staggered magnetization versus temperature in 
YBa2Cu30y for several values of doping. Values of doping 
X = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 are shown near the correspond- 
ing curve. These curves have been obtained without account 
of the lifetime of the hole, F = 0. 

presented in Fig. [5] are in a qualitative and to some ex- 
tent quantitative agreement with data from Refi^ shown 
in Fig. [Tfbottom) and in Fig. |9l[bottom). The theory 
clearly demonstrates that while reductions of the zero 
temperature staggered magnetization at x < xq is pretty 
small. Fig. [TJtop), the reduction of the Neel tempera- 
ture with doping is dramatic. There are two reasons for 
the reduction, (i) Thermal excitation of the precursor to 
the hourglass, the lower incoherent part of the magnetic 
spectrum shown in Fig. [5] (ii) Thermal population of 
the 7 = — 1 band. Fig. |8] indicates also some negative 
bending of /i(r) curves in a qualitative agreement with 
data presented in Fig. [SJbottom). 

Theoretical curves plotted in Fig. |8] demonstrate even 
too steep decrease of the Neel temperature with doping 
compared to experimental data shown in Fig.[IJbottom). 
For example at a; = 0.03 the theoretical Neel temperature 
is 140K, Fig. [H while experimentally it is about 300K, 
Fig. [IJbottom). To fix this problem one has to realize 
that the above consideration of ferniions disregards an 
important physical effect, the finite lifetime (scattering 
time) of a fermion at a nonzero temperature. This is 
the same effect that leads to the temperature dependent 
resistivity. To understand importance of this effect prior 
to calculations one has to recall, see previous subsection. 



that the Neel temperature is formed at very large in-plane 
distances up to 100 lattice spacing. This corresponds to 
q ^ 0.01 in Eqs. (pi)) . (P5|) . The fermionic polarization 
operator at so small q (so large distances) must depend 
on the fermion mean free path. This explains crucial 
importance of the fermion lifetime. To account for the 
lifetime effect Eq. ([^S]) has to be modified in the following 
way 

^o(^c,.) = 2.AcW E E-4n?r^^''^ 

Here F is the broadening due to scattering. Let us take 
the usual width characteristic for the two-dimensional 
Fermi liquid^, 

r = A— , (29) 

where ep — /3p|,/2 — i:(3x is the Fermi energy. I disre- 
gard the logarithmic T-dependence of the coefficient A, 
the dependence is beyond accuracy of the calculation. 
The coefficient A will be used as a fitting parameter. 
Note, that generally the width F depends on both tem- 
perature T and Matsubara frequency . The dominating 
contribution to Ea. p4)) comes from the zero Matsubara 
frequency. Therefore, the width F is important in the 
zero frequency term, s = 0, and it is completely negligible 
in s 7^ terms. Hence, the width ((29)) corresponds to the 
zero Matsubara frequency. Numerical evaluation of the 
polarization operator (|28p is not more difficult than eval- 
uation of (pSj) . Solution of RG equations gives the stag- 
gered magnetization ^.{T, x) with account of the fermion 
lifetime. The best fit of the experimental dependence of 
the Neel temperature on doping is achieved at 

Aw 0.7. (30) 

The calculated Neel temperature versus doping is shown 
by large ired dots in Fig. [TJbottom). The calculated 
staggered magnetization versus temperature is plotted 
in FiglHKtop) for several values of doping. Experimental 
curves from Ref. |^ are presented in Fig. IHKbottom) . Over- 
all agreement between theory and experiment is very 
good. 

It is worth stressing again that the calculation of the 
temperature dependence of the magnetization in the lay- 
ered system is less reliable than calculations of zero tem- 
perature properties in section HI. The complexity of the 
finite temperature case is due to the very large span of 
spacial scales involving in the problem with the largest 
scale about 100 lattice spacing. Only leading effects 
have been taken into account in the present calculation. 
Glearly there are subleading effects that also influence the 
magnetization. For example, usual disorder (impurities) 
must influence fermion dynamics at the scale ~ 100 lat- 
tice spacing and hence influence magnetization. In view 
of this comment the agreement between the theory and 
experimental data, see Fig. [IJbottom) and Fig. [SI is re- 
markable. Most importantly, the theory explains why the 
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80 160 240 320 400 480 

FIG. 9: Staggered magnetization versus temperature 
in YBa2Cu30j; for several values of doping x. The- 
oretical curves with account of the hole lifetime for 
X = 0,0.1,0.2,0.3,0.4,0.5,0.6 are shown in the top fig- 
ure. Experimental curves from Ref. are presented 
in the bottom figure, the doping levels are x = 
0.001, 0.002, 0.005, 0.02, 0.029, 0.036, 0.039, 0.049, 0.061, 0.065. 



Neel temperature drops down dramatically with doping, 
while the zero temperature magnetization is almost dop- 
ing independent. This "contradictory" behaviour is due 
to the band splitting, and due to different fillings of the 
split bands. The different filling is a fingerprint of small 
hole pockets. At zero temperature only the lower band 
is occupied while temperature populates the upper band 
as well. The "contradictory" behaviour is closely related 
to the Lifshitz point at a; = ccg and to the development of 
the spin spiral aX x > xq when both bands are occupied 
at zero temperature. 



V. CDW INDUCED BY OXYGEN CHAINS, 
SMALL HOLE POCKETS, MECHANISM FOR 
BAND SPLITTING 

The key point of the YBCO phenomenology suggested 
in Ref. |lj and applied in the present work is the splitting 
of hole bands. There are other key points like small hole 
pockets, spin spirals, etc. However these other points are 
not specific to YBCO, they are generic for aU cuprates. 
The band splitting is specific to YBCO. The paper— sug- 
gested that the hole band splitting in YBCO was due 
to the hole hopping between layers inside the bilayer, 
the bonding-antibonding splitting. However, our recent 
analysis^ indicates that antiferromagnetic correlations 
between the layers forbid the bonding-antibonding split- 
ting. If the hole hopping matrix element between the 
layers is t± then the band splitting in the case of AF 
correlations between the layers is oc t^icoskx + cosfcy), 
see Ref. [U. The splitting is zero at the nodal points 
{kx,ky) = (±7r/2, ±7r/2) contrary to the assumption ([T]). 




FIG. 10: Single hole dispersion in the AF background. 



Thus, contrary to the assumption in Ref. uA the inter- 
layer hopping cannot contribute to Aq. Another mecha- 
nism for splitting is necessary. In this section I argue that 
the splitting is due to oxygen chains. Let us first consider 
the case y = 6.5 when every second chain is fully filled. 
In this case chains produce the effective pscudopotcntial 
for in-plane holes 



VciX) = -?;o cos(7rA:) 



(31) 



where Vq is the amplitude of the potential and X is the 
direction perpendicular to chains (I denote the distance 
by capital X to make it different from the doping x). 
The hole dispersion in the antiferromagnetic background 
is shown in Fig. 1101 The dispersion consists of two full 
pockets, the pocket a and the pocket b. There is the 
perfect nesting condition between the pockets and the 
chain potential (|3T|) . therefore the two split bands are 
generated 



2 , 
2 , 

-'■'0 + ' V- 



|q)p + 
^/2 



(X cos 



(|a) e-' 



oc sin e^P ' 



(32) 



Here p is the momentum with respect to the center of 
the pocket. Eqs. ([5^ represent exactly the 7 = ±1 
band splitting adopted in ([7|). Due to the exact nesting 
of small hole pockets a tiny pseudopotential vq ~ 30meV 
is sufficient to generate the splitting Aq = 2vq « QbmeV 
that follows from the magnetic analysis, see Eq.®. 

There might be an impression that the splitting ((32)) 
is not quite equivalent to the effective action ([7]). For 
example, the question arises why there is no a spin-wave 
vertex that transfers ip^ to tp^ ? The vertex carries a large 
momentum tt, therefore the vertex vanishes after integra- 
tion over X . In other words, soft magnons included in the 
effective action ([7]) cannot induce a transition with the 
large momentum transfer. A careful analysis shows that 
the splitting ([5^ with account of two layers is completely 
equivalent to ([7]). 

According to ([32)1 the wave function of the lower 7 = 
-1-1 band is nonzero at X=0,2,4,... while the wave func- 
tion of the upper 7 = —1 band is nonzero at X=l,3,5,... 
Due to the splitting the bands are differently populated 
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and this results in the in-plane hole density wave with 
period of the two lattice spacingJ^ii^ Let us calculate 
the amplitude of the CDW. The oxygen content y = 6.5 
corresponds to doping x « 0.1, see Eq.(IT|). This doping 
is within the spin-spiral phase, therefore to calculate fill- 
ings of bands one has to account the spin spiral as it is 
shown in Fig. Upright). In the lower 7 = -|-1 band both 
pseudospin projections are populated, while in the upper 
7 = — 1 band only one pseudispin projection is populated. 
Populations of different subbands have been calculated in 
the analysis of magnetic quantum oscillations, see Eqs. 
(4) in Ref. [H. From these equations one concludes that 
populations of the upper and lower bands (7 = =pl) are 



A 



2(3 -2A) 

~ X — x^i . 



{x - Zo) 



(33) 



Naturally, the population of the upper band vanishes at 
x = .To, this is the Lifshitz point. For a; = 0.1, Xq = 0.06, 
and A — 1.23, one finds x-i = 0.03 and x+i = 0.07. 
Hence, the hole density per cite at every even value of 
X is 2a;+i = 0.14, and at every odd value of X it is 
2x_i — 0.06. However, this is not the amplitude of the 
CDW yet. 

All equations in the present paper arc written in terms 
of holes dressed by magnetic quantum fluctuations (mag- 
netic polarons). Hence x±i are densities of the dressed 
holes. The dressed hole has a finite size, therefore, the 
charge density modulation is smaller than that naively 
given by x±i. It is known that the quasiparticle residue 
of the dressed hole is about Z « 0.4, see e.g. Ref. HI. 
This means that with the probability Z ^ OA the hole 
resides at the same site as the quasihole and with the 
probability (1 — Z)/A « 0.15 the hole resides at each 
of the four nearest Cu sites. Therefore, the real charge 
densities per site are 



P+i = 2 
P-i = 2 



l-Z\ 1~Z 
Z + 2 — - — x+i + 2 — - — x^i 
4 / 4 

1-Z\ 1-Z 
Z + 2 — : — I x-i + 2 — - — x+i 



0.12 

0.08 

(34) 



This gives the amplitude of the CDW. The estimate 
of the amplitude is based purely on magnetic data, it 
depends mainly on the position of the Lifshitz point, 
xo = 0.06. 

NQR was not used in the estimate (p4)). Neverthe- 
less the estimate is pretty much consistent with NQR 
datapi^ presented in Fig. |31 The NQR frequency shift 
with respect to the frequency of the undoped sample, 
I'o = 23.3MHzji^ is proportional to the local hole con- 
centratior>ii 



VQ = 23.3MHz + Bp 



(35) 



The higher frequency NQR line at 1/ = 6.56 is ~ 
30.3MHz, see Fig. ||l According to Eq.dM]) the hue 



corresponds to p « 0.12. Hence, the constant B in 
Eq.dSSl) is B « bSMHz/hole. Interestingly, the value 
of B is significantly larger than that in La2-xSra;Cu04, 
B « 20MHz/hole^ and in HgBa2Cu04+5, B « 
30MHz/hole^ Assuming that optimal doping corre- 
sponds to p w 0.14 and using Eq. (P5|) one finds the opti- 
mal doping NQR frequency ly = 31.AMHz. This value is 
pretty close to i^opt ~ 31.6M_ffz that follows from Fig. [3] 
at ?/ w 7. According to (|34)) the lower frequency NQR 
line at y sa 6.5 corresponds to p « 0.08. Substituting this 
value in Eg. ([55]) one finds the frequency v = 27.9MHz. 
Again, this value is pretty close to the lower frequency 
line i/i sa 27.8MHz that is shown in Fig. [3]at y = 6, 56. 
Thus, the amplitude of the CDW determined from the 
position of the Lifshitz point is fully consistent with the 
NQR data. 

The simple potential (pi]) is literally applicable only 
to y = 6.5. Obviously, there is no any modulation at 
y ~ 7 as the rightmost Fig. [3] indicates. Away from y = 
6.5 more complex oxygen superstructures can appeari^ 
Assumption important for the present work is that at 
< X < 0.1 (6.20 < y < 6.5) the superstructure ((3T|) 
is dominating. NQR datapi^ for y = 6.4 and y = 6.45 
presented in Fig. |3] support this assumption: there are 
only two NQR lines that are only slightly broader than 
the lines at y = 6.5. 



VI. CONCLUSIONS 

Small hole pockets and the associated spin spiral state 
are generic properties of all cuprates at low doping. The 
key point of the YBCO phenomenology additional to the 
generic properties is splitting of the hole pockets into the 
lower band and the upper band. This splitting together 
with splitting of magnon to the acoustic and the optic 
mode provides stability of the collinear anifcrromagnetic 
order at doping below the Lifshitz point at x ~ 0.06. 
At doping below the Lifshitz point only the lower band 
is populated. At higher doping the upper band starts 
to populate and simultaneously the spin spiral starts to 
develop. 

At doping below the Lifshitz point the doping induced 
spin quantum fluctuations are pretty weak. This ex- 
plains why the zero temperature staggered magnetization 
is close to 0.6ps, the maximum value allowed by quan- 
tum fluctuations of localized spins. The developed theory 
perfectly reproduces the weak decrease of the staggered 
magnetization with doping observed experimentally. 

While the zero temperature staggered magnetization 
is almost doping independent, the Neel temperature de- 
cays very quickly from Tn = A20K at x = to practically 
zero at X sa 0.06. This quick decay is a consequence of 
the closeness to the Lifshitz point. Again, the theory 
reproduces very well the doping dependence of the Neel 
temperature as well as the observed temperature depen- 
dence of the staggered magnetization at a given doping. 

The band splitting (the hole pocket splitting) is in- 
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duced by the modulation of oxygen chains. The main 
period of the modulation is two lattice spacing. Because 
of the perfect nesting between the small hole pockets 
and the period of the modulation, a small pseudopoten- 
tial caused by the chains is sufficient to induce the band 
splitting about 60meV. The splitting causes the in-plane 
charge density wave with a significant amplitude depen- 
dent on doping. 
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